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Abstract 

We study anomalous diffusion for one-dimensional systems described by a generalized 
Langevin equation. We show that superdiffusion can be classified in normal superdiffusion 
and fast superdiffusion. For fast superdiffusion we prove that the Fluctuation-Dissipation 
Theorem does not hold, which induces an effective temperature in the system. This effective 
temperature is a signature of metastability found in many complex system such as Spin-Glass 
and granular material. 

1 Introduction 

Since its formulation, the Fluctuation-Dissipation Theorem (FDT) has played a central rolepJE] 
in non-equilibrium statistical mechanics in the linear response regime (LNESM). It reaches such 
an importance that a full formulation of LNESM is given j2] based on it. In the last 30 years, 
fundamental concepts and methods have been developed P] El IH HI ED and a large number 
of connections have been established. A necessary requirement for the FDT is that the time- 
dependent dynamical variables are well defined at equilibrium. The presence of nonlinear effects 
or far from equilibrium dynamics may lead to situations where the FDT does not hold, the aging 
process in spin-glass systems being a good example [71 iRl Ifll fTol fTTl fT2l fl3j . 

In this letter, we present a straightforward proof of the inconsistency of the FDT for a 
certain class of superdiffusive processes described by a generalized Langevin equation (GLE). 
The surprising result here is the fact that the violation happens in the linear regime, i.e. where 
we expect the validity of the LNESM. The use of the FDT allows us to classify two classes of 
superdiffusion. The first class, which we shall call normal superdiffusion, does obey the FDT; 
the second class, which we shall call fast superdiffusion, does not obey the FDT. The proof is 
simple and we discuss as well how the diffusive process leads to an equilibrium. 

Diffusion is one of the simplest processes by which a system reaches equilibrium. For normal 
diffusion, the process is so well known that it may be described by an equilibrium type distri- 
bution for the velocity and position of a particle. However, the strange kinetics of anomalous 
diffusion, intensively investigated in the last years |14U15lH7illl7lH%lllfl| . shows surprising results. 
Consequently, studying anomalous diffusion seems to be a good way to obtain the conditions of 
validity for the FDT. 
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2 Violation of the FDT 

We shall start writing the GLE for an operator A in the form pfl Ej 



dA{t) = - f T(t - t')A(t')dt' + F(t), (1) 



dt 

where F(t) is a stochastic noise subject to the conditions (F(t)) = 0, (F(t)A(0)) = and 

C F (t) =< F(t)F(p) >=< A 2 > eq r(t). (2) 

Here CV(i) is the correlation function for F(t) and the brackets <> indicate ensemble average. 
Eq. Q is the famous Kubo FDT and it is quite general^]. In principle, the presence of the 
kernel T(t) allows us to study a large number of correlated processes. The main quantity is the 
correlation function CU(£) =< A(t)A(0) >, from which we can describe most of the process of 
interest. 

We may naively expect that, by Eq. Q and Eq. (2J, a system will be driven to an equilibrium 
, i.e. 

lim < A 2 (t) >=< A 2 > e „ . (3) 

t— >oo 

We shall see however that this is not always the case for superdiffusive dynamics. Let us define 
the variable 

y(t)= f A(t')dt', (4) 
Jo 

with asymptotic behavior 

lim < y 2 (t) >~t a . (5) 

t — >oo 

For normal diffusion a = 1, we have subdiffusion for a < 1 and superdiffusion for a > 1. 
Notice that if A(t) is the momentum of a particle with unit mass, y(t) is its position. Using the 
Kubo's definition of the diffusion constant, Morgado et al|32| obtained a general classification 
for anomalous diffusion; i.e. using 

D = lim C A (z) = lim < A \ >eq , (6) 
«-o Ay ' z^o Z + T(z) V ' 

where T(z) is the Laplace transform of T(t) they obtained: For a finite value of T(0) ^ 0, normal 
diffusion; for T(0) = 0, superdiffusion; and for T(0) = oo, subdiffusion. 
Now the Laplace transform of Eq. Q suggests a solution of the form 

A(t) = f R(t- t')F(t')dt', (7) 



where we have set A(0) = 0, and 



- — ipr- (8) 

z + T(z) 



Squaring Eq. (JJ and taking the ensemble average we obtain 
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< A 2 (t) > = 2 / R{t') ( C F (t' -t")R{t")dt"dt'. (9) 
Jo Jo 

At this point, it is quite usual to perform numerical calculation [Uj. However, we shall show 
here that very important results can be obtained analytically. From Eq. <|HJl, we can get a 
self-consistent equation for R{t) as 

= - f T{t-t')R{t')dt'. (10) 

dt In 



Notice from Eq. © and Eq. © that R(t) is the normalized correlation function, i.e. R(t) = 
Ca(£)/Ca(0). By using the FDT Eq. (|2j and Eq.^Ol we can exactly integrate Eq. © and 
obtain 



where 



< A 2 (t) >=< A 2 > eq A(t), (11) 
A(i) = \-R 2 {t). (12) 



Equation. is satisfied if and only if 

lim X(t) = A* = 1, (13) 

t — >oo 

or equivalently 

lim R(t) = lim zR{z) = lim % = 0. (14) 

t->oo z^o z^o z + r(z) 

Where in the first step we used the final value theorem[23|. Again we shall comment here that 
if Eq. l[T4Tl does not hold, its equivalently Eq. does not hold as well, and it violates the FDT 
Eq. (2J. Recently Lee[2j, using his recurrence relation formalism^, obtained a condition similar 
to Eq. ltl4Tl for the ergodicity hypothesis to work. Fe. if the system keeps some memory after an 
infinity time, the ergodicity does not hold. We show here that it is the same condition for the 
validity of the FDT. 

Equation l(T4Tl is satisfied for normal and subdiffusion. For superdiffusion \im z ^ T(z) = 0, 
and Eq. ltl4Tl becomes 

lim J R(t) = (l + limffi)- 1 , (15) 

t-+oo z^o az 

There are two distinct limits for Eq. 11511 . which define two classes of superdiffusion. For the first 

class, lim z ^ dF g^ = 00 an ^ tne system obeys the FDT. The second class has lim z ^ ^= oo 
and it does violate the FDT. The first class we shall call normal superdiffusion (NSD) and the 
second class fast superdiffusion (FSD). 

Consider now the asymptotic behavior for T(z)as 

f (z ->• 0) = az". (16) 

Is easy to see that for v < we have subdiffusion, for v = normal diffusion, and for v > 0, 
superdiffusion. From the above equations we have for < v < 1 NSD and, finally, for v > 1 we 
have FSD. There is an obvious connection between v and a, defined in Eq. 0, that classify the 
diffusion. Morgado et al [IJI] shown that v = a — 1 and consequently the FSD starts at a > 2, 
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i.e., the FDT does not work for the ballistic motion and beyond. Better, the formalism of GLE 
+ FDT works for < a < 2. We shall discuss later the lower limit a — 0. 



3 Ballistic Motion 

The ballistic motion is at the limit of the validity of the FDT, to go beyond it is too dangerous so 
we shall keep working within a — 2. Before we give some example let us make a very important 
association. The force F(t) in Eq. Q can be obtained from a thermal bath composed of harmonic 
oscillations, consequently accordingly to Eq. <|2jl the memory can be put as 

T(t) = J p(ui)cos(u)t)du), (17) 

where p{u) is the noise density of states. The same argument used before for the Laplace 
transform can be used for for the Fourier transform, with the simplifying consequence T(oj) = 
p{oj). This is a great advantage, since noise density of states exist not only for systems governed 
by GLE, but for most of the physical systems. Consequently the MOBH conjecture |33| reads: // 
a disordered unidimensional system has p(u) ~ u> u , as ui — > than the diffusion exponent is 



a = v + l. (18) 

This conjecture has been observed for the quantum disordered Heisenberg ferromagnetic 
chain[2I], and is under discussion for energy propagation on the harmonic disordered chain[22j. 
We can choose now the density of state to produce the diffusion we want. Consider now the 
noise density of states 

, . f constant u\ < oj < u>2 
p{uj) = { otherwise ' (19) 

for u>i — we have the Debye density of states for a thermal noise made out of acoustic phonons. 
Thus, for u>i — we have normal diffusion and for any wi ^ we have superdiffusion. This 
density yields 



r(t) = 



sm(uj2t) sin(wit) 



t t 



(20) 



The Laplace transform of Eq. l|20j) gives as z — > T(z)~z, consequently v = 1 and a = 2, which 
is the ballistic limit. If we let (3 — c^ 2 /2 we get A* as 

( 2ljJ-\LjJ2 \ 2 , 

A = 1 - I — . 21 

Values of A*^l show the inconsistency of the FDT because we start supposing the existence of an 
equilibrium value < A 2 > eq and, after an infinite time, we end up with < A 2 > eq A*. Equation 
1)2111 has a parameter control u>i, which measures the "hole" in the density of states, and how far 
we are from the result predicted by the FDT. 

Now we select A(t) = v(t), the particle's velocity, so that < v 2 (t) >=< v 2 > eq X(t). We 
simulate the GLE for a set of 10, 000 particles starting at rest at the origin and using the 
memory in Eq. 12011 with u>2 — 0.5 and different values of u>\. The results of these simulations 
are shown in Fig. 1, where we plot < v 2 (t) >. We used the normalization < v 2 > eq = 1, so that 
< v 2 (t) >= X(t). Notice that X(t) does not reach a stationary value, rather it oscillates around 
a final average value A s . This value of A s should be compared with A* obtained from Eq. l(2lll . 
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Figure 1: Normalized mean square velocity as a function of time for the memory given by Ea. l|20ll . 
Here (3 = UJ2/2 and u>2 = 0.5. Each curve corresponds to a different value of u>\. a) uj\ = 0; 
b) uj\ — 0.25; c) u)\ = 0.45. The horizontal lines correspond to the final average value A s . In 
agreement with the theoretical prediction, As decreases as u\ grows. 

Figure 2: A* as a function of the parameter W\. Each dot corresponds to a value of A s obtained 
from simulations like those described in Fig. 1. The line corresponds to the theoretical prediction 
given by Eo. l(2ll . 

In Fig. 2 we plot A* as a function of u\ as in Eq. 1I2TI) with fixed u 2 — 0.5. We also plot 
the final average values \ s obtained from simulations for different values of u>\ . Notice that as 
u>i increases A* decreases as expected. The agreement between simulations and Eq. 12 "H shows 
that we can predict the average value X s , even when the FDT does not work. 

Now we can define A*= T*/T, where T* is an effective temperature for the system. Effective 
temperatures different of the expected temperature T, or A*^ 1, are found in spin glasses where 
the FDT does not work[H| El CHI EH- The first observation of such phenomena was reported 
Kauzmann0. He noticed that if the entropy of a supercooled liquid is extrapolated below 
the glass temperature T g , it becomes equal to the crystal temperature T c > 0, and in some 
cases even T c < 0. To avoid this paradox he suggested the existence in the supercooled liquid 
phase of effective spinodal temperature T sp <T C . In a recent work Rubi et al|3} investigate the 
violation of the FDT, using a Fokker-Planck approach. They found temperatures T* which 
are greater and small than T. Ricci-Tersenghi et al[10j and Cavagna et al llj performed single- 
spin-flip Monte Carlo simulations in square lattices with frustration and they obtained effective 
temperatures T* ^ T. Methods for measuring those effective temperatures are discussed as well 
in the literature^]. 

An oscillatory behavior similar to that found in Fig. 1 was observed by Srokowski fTTl Il8| . 
in his simulations using GLE. However, the kind of motion he studied is a subdiffusive motion. 
Using his memory we get^HI = j3[l — exp(ez)/(ez) + Ei(ez)], here f3 is a constant, e a small 
number, and E\ is the exponential integral. This result gives lim z ^ T(z) ~ z^ 1 , from the above 
arguments we get v = — 1 and a = 0. We believe we can explain this strange results. Let us 
consider a constant memory of the form T(t) = K, if A(t) in Eq. Q is the velocity, than the 
memory term yields —Ky, where y is the position. Consequently the particle is bound to the 
origin by a harmonic spring and has no diffusive behavior. Indeed the Laplace transform gives 
r(z) = K/z, with v = — 1 and a = 0, the same as the result we obtain from Srokowski memory. 
We shall notice that a — does not meaning that the motion is a harmonic type motion, rather 
it means it belongs to same class, the diffusion behavior of Srokowski is not a power low as in 
Eq. Q probably it is slower than a power low, such as, for example, a logarithmic behavior. 

For a > 2, the FSD cannot be described by the methods we used here. Once the FDT does 
not work, the GLE and the FDT together predict strange results such as a null dispersion for the 
dynamical variable, i.e. < A 2 (t — > 00) >= 0. Moreover, the exponent a can be put as a = 2/ Dp, 
where Dp is the fractal dimension [22]. Consequently a > 2 leads to Dp < 1, which is not a full 
curve, but a set of points such as the Cantor set, and cannot represent a classical trajectory. 
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4 Conclusion 



We discussed the stationary behavior for the mean square value of a dynamical variable A(t) 
and as well the mean square displacement of the quantity' y(t) — J Q A(t')dt' . The assymptotic 
behavior of < y 2 (t) >~ t a as t — > oo, can be explained for < a < 2. We show the the 
superdiffusive motion must be classified in normal superdiffusive (NSD), for 1 < a < 2, and 
fast superdiffusive (FSD), for a > 2. The FSD motion shows an inconsistency between the 
GLE and the FDT. This kind of superdiffusion with a > 2 is common in hydrodynamical 
processes. It is not surprising that these processes will be far from equilibrium and violate the 
FDT. We pointed out here how it happens and precisely where the FDT breaks down. As we have 
already mentioned, spin glasses seem to be a rich field for studying these phenomena. Indeed 
experimental an d theoretical works E3 EH have been reported in this area, confirming 
the violation of the FDT. As well, the effective temperature found in noncrystaline material 
is connected here with the FSD and the violation of the FDT. It would be very helpful if the 
exponent a for those diffusive processes could be measured. Another related phenomena are 
anomalous reaction rate[24J and chaos synchronization 2Sj, which we expect to discuss soon. 
Although anomalous diffusion remains as a surprising phenomena, we hope that this work will 
help in the centennial effort to understand diffusion and the relation between fluctuation and 
dissipation. A generalization of the FDT to include the FSD is necessary, what will require a 
deeper understanding of systems far from equilibrium. 

References 

[1] R. Kubo, Rep. Prog. Phys. 29, 255 (1966). 

[2] R. Kubo, M. Toda, N. Hashitsume, Statistical Physics II, Springer (1991). 
[3] H. Mori, Prog. Theor. Phys., 33, 423 (1965). 

[4] D. J. Evans and G. P. Morriss, Statistical Mechanics of Nonequilibrium Liquids, Academic 
Press (1990). 

[5] M. H. Lee, J. Math. Phys., 24,2512 (1983)-Phys. Rev. Lett, 85, 2422 (2000) 
[6] M. H. Lee, Phys. Rev. Lett., 87, 601 (2001). 
[7] G. Parisi, Phys. Rev. Lett., 79, 3660 (1997). 
[8] W. Kauzmann, Chem. Rev. 43, 219 (1948) 

[9] I. Santamaria-Holek, D. Requera and J. M. Rubi Phys. Rev. E, 63, 051106 (2001) - A. 
Perez-Madrid, D. Requera and J. M. Rubi, cond/mat 0210089. 

[10] F. Ricci-Tersenghi, D. A. Stariolo and J. J. Arenzon, Phys. Rev. Lett, 84, 4473 (2000). 

[11] A. Cavagna, I. Giardine and T. Grigera, Europhys. Lett., 61, 74 (2003). 

[12] R. Exatier and Peliti, Eur. Phys. J. B 16, 119 (2000). 

[13] T. S. Grigera, N. E. Israeloff, Phys. Rev. Lett. 83, 5038 (1999). 

[14] J. P. Bouchaud and A. Georges, Physics Reports, 195, 127 (1990). 

[15] H. Scher, M. F. Shlesinger and J. T. Bendler, Physics Today, January, 44, 26 (1991). 



6 



[16] M. F. Shlesinger, G. M. Zaslavsky and J. Klafter, Nature, 363, 31 (1993). 

[17] T. Srokowski, Phys. Rev. Lett., 85, 2232 (2000)-P%s. Rev. E, 64, 31102 (2001). 

[18] F. A. Oliveira, R. Morgado, C. Dias, G. G. Batrouni andA. Hansen, Phys. Rev. Lett., 86, 
5839 (2001). 

[19] R. Morgado, F. A. Oliveira, G. G. Batrouni and A. Hansen, Phys. Rev. Lett, 89, 100601 
(2002). 

[20] M. R. Spiegel, Theory and Problems of Laplace Transforms, Schaum Publishing Co. (1965), 
page 20 . 

[21] F. A. B. Moura, M. D. Coutinho, E. P. Raposo, Phys. Rev. B, 65 104416 (2002). 

[22] F. A. B. Moura, M. D. Coutinho, E. P. Raposo, cond-mat/0209304. 

[23] D. Stauffer and H. E. Stanley, From Newton to Mandelbrot, Springer (1995). 

[24] F. A. Oliveira, Physica A 257, 128 (1998); A. M. Maroja, F. A. Oliveira, M. C. Ciesla and 
L . Longa, Phys. Rev. E, 63, 061801 (2001) 

[25] M. C. Ciesla, S. P. Dias, L. Longa and F. A. Oliveira Phys. Rev. E., 63, 5202 (2001); L. 
Longa, E. M. Curado and F.A. Oliveira, , Phys. Rev. E, 54, 2201 (1996). 



7 



0.8 - 



^ 0.6 - 




0.4 - 



0.2 - 



CO 



